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D ! Abstract 

. We study the evolution of a two dimensional minisuperspace eosmological model in classical 

^>0 ' and quantum levels by the Noether symmetry approach. The phase space variables turn out to cor- 

respond to the scale factor of a Friedmann-Robertson- Walker (FRW) model and a scalar field with 
which the action of the model is augmented. It is shown that the minisuperspace of such a model 
O , is a two dimensional manifold with vanishing Ricci scalar. We present a coordinate transformation 

which cast the corresponding minisuper metric to a Minkowskian or Euclidean one according to 



^—^ . the choices of an ordinary or phantom model for the scalar field. Then, the Noether symmetry of 

such a eosmological model is investigated by utilizing the behavior of the corresponding Lagrangian 
under the infinitesimal generators of the desired symmetry. We explicitly calculate the form of 
' the scalar field potential functions for which such symmetries exist. For these potential functions, 

, the exact classical and quantum solutions in the cases where the scalar field is an ordinary or a 

ly-^ ' phantom one, are presented and compared. 

m 
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1 Introduction 

^ \ Classical, semiclassical and quantum scalar fields have played a central role in conceptual discussion 

' of unified theories of interactions and also in all branches of the modern eosmological theories. From 

a eosmological point of view, there is a renewed interest in the scalar-tensor models in which a non- 
minimal coupling appears between the geometry of space-time and a scalar field [1]. This is because a 
number of different scenarios in cosmology such as spatially flat and accelerated expanding universe at 
the present time [2], inflation [3], dark matter and dark energy [4], and a rich variety of behaviors can 
be accommodated phenomenologically by scalar flelds. Traditionally, eosmological models of inflation 
use a single scalar field with a canonical kinetic term of the form l/2g'^'^d^(j)du(j) with some particular 
self-interaction potential V{(l)) like l/2m(/)^ or A(/>'^, etc. Such a scalar field is often known as minimally 
coupled to the geometry. However, there are also scalar fields in what is qualified to be called the 
scalar-tensor theory, which are not simply added to the tensor gravitational field, but comes into play 
through the non-minimal coupling term [5]. On the other hand, the current eosmological observations 
allow the possibility of existence of a cosmic fluid with equation of state parameter smaller than — 1 
[6]. One of the simplest explanations of such an equation of state is considering a scalar fleld with 
negative kinetic energy which is called a phantom field [7]. The classical and quantum eosmological 
dynamics of the phantom scalar fields has been studied in a number of works, see for example [8]. 

In this paper we shall study the classical and quantum dynamics of a flat FRW model with a scalar 
field as its source by the Noether symmetry approach. We set up an effective Lagrangian in which the 
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scale factor a and scalar field (f) play the role of independent dynamical variables. This Lagrangian is 
constructed in such a way that its variation with respect to a and (f) yields the appropriate equations 
of motion. The form of the potential function of the scalar field is then found by demanding that 
the Lagrangian admits the desired Noether symmetry [9]. By the Nocthcr symmetry of a given 
minisuperspace cosmological model we mean that there exists a vector field X, as the infinitesimal 
generator of the symmetry, on the tangent space of the configuration space such that the Lie derivative 
of the Lagrangian with respect to this vector field vanishes. In [10] the applications of the Noether 
symmetry in various cosmological models are studied. Here, before applying the Noether symmetry 
condition on the model under consideration, we shall see that the corresponding minisuperspace of our 
model is a two dimensional Riemannian manifold with zero Ricci scalar. Therefore, it is possible to 
find a new set of coordinates in terms of which the minisuper metric takes the form of a Minkowskian 
(when the scalar field is an ordinary quintessence field) or Euclidean (when the scalar field is a 
phantom field) space. The desired coordinate transformation brings the Lagrangian into a canonical 
form where its kinetic and part of its potential terms are like the ones of a couple of harmonic 
oscillators. We shall see by demanding the Noether symmetry as a feature of this Lagrangian, we 
can obtain the explicit form of the potential function. Since the existence of a symmetry results in a 
constant of motion, we can integrate the field equations which would then lead to expansion law of 
the universe. 

The structure of the article is as follows. In section 2, we briefly present the basic elements of 
the issue of Noether symmetry. In section 3, we first introduce the scalar field FRW cosmology 
and write its Lagrangian in terms of the minisuperspace variables. This section is then divided into 
two subsections each of which will deal with one kind (quintessence or phantom) of scalar field. In 
these two subsections, after applying the Noether symmetry condition on the model, we provide 
the analytical solutions for the corresponding Noether symmetric classical and quantum cosmologies. 
Finally, the conclusions are summarized in section 4. 

2 The Noether symmetry 

In this section we briefly review the issue of the Noether symmetry approach used in the present work 
which originally appeared in [11]. To do this, we consider a dynamical system with finite degrees of 
freedom moving in a Riemannian space with the metric components Qab- The evolution of such a 
system can be produced by the following action 



where q are the dynamical variables representing the configuration space of the system (the indices A, 
B, ... run over the dimension of this space), U (q) is the potential function and r is an affine parameter 
along the evolution path of the system. In time parameterized theories, for instance general relativity, 
in which their action is invariant under time reparameterization, the affine parameter r can be linked 
to a time parameter t by a lapse function N(t) through Ndt = dr. Therefore, for such a dynamical 
systems the action (1) can be written as 



where now an over dot indicates derivation with respect to the time parameter t and £(q, q) is 
the Lagrangian function of the system. Variation of the above action with respect to yields the 
equations of motion as 




(1) 



5 = ^ dtCiq^, q^) = dtN y^g^Bq^q^ - W(q)] , 



(2) 




(3) 
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where = -g'^^dsU is a force term and T^^ = ig^^ {ObGdc + OcQbd - doGBc) is the Christof- 
fel symbols associated with the metric Qab- In this formaHsm N is not a dynamical variable in the 
sense that variation of the action (2) with respect to it yields 



2iV2 



^AB^V+W(q) = 0, 



(4) 



which is nothing but the Hamiltonian constraint, a constraint that should hold during the evolution 
of the system and is indeed a reflect of the fact that the underlying theory is a time parameterized 
theory. Basically, the equations of motion (3) are the geodesic equations of a "point particle" (=the 
system) moving in a Riemannian space under the act of the potential field U{c\). Quantization of 
such a dynamical system may be achieved by the method of canonical quantization. For this purpose 
we introduce the momenta conjugate to the variables q"^ by 



Pa = 



dc 



Gab 



iV' 



leading to the following Hamiltonian function 



H = Paq^ -C = N 



-g^^'PAPB+Uici) 



= Nn. 



(5) 



(6) 



Because of the constraint equation (4), the above Hamiltonian seems to be identically equal to 
zero. Therefore, under canonical quantization this Hamiltonian yields the Wheeler-DeWitt (WDW) 
equation 7{*(q) = 0, where ^'(q) is the wave function of the quantized system and T-L should be 
written in a suitable operator form. If one makes a natural choice of factor ordering, the WDW 
equation may be written as [12] 



H*(q) = 



1 



V2+W(q) 



*(q) = 0, 



(7) 



where = -^^^Oa {V—QQ^^dB^ is the Laplacian operator in the space with metric Qab- 

Now, let us deal with the idea of the Noether symmetry in a given dynamical system like that one 
presented above. Following [11], we define the Noether symmetry as a vector field X on the tangent 
space TQ = (q, q) of the configuration space through 



dt dq^ ' 



(8) 



where a'^(q) are unknown functions on configuration space. The Noether symmetry then implies 
that the Lie derivative of the Lagrangian with respect to this vector field vanishes, that is, LxJC = 0, 
which leads 



a 



dC da^{q) dC 



+ 



dt dq 



= 0. 



(9) 



Noether symmetry approach is a powerful tool in finding the solutions to a given Lagrangian, including 
the one presented above. In this approach, one is concerned with finding the cyclic variables related 
to conserved quantities and consequently reducing the dynamics of the system to a manageable one. 

The existence of Noether symmetry means that phase flux is conserved along the vector field X and 



thus a constant of motion exists. 

Euler-Lagrangc equations = 
motion are found as 



Indeed, noting that Pa 
from (9) we get 4 



dC 



dC 



and also taking into account the 



a' 



\<i)Pa 



Q = a^(q)PA. 



0. Thus the constants of 



(10) 
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In order to obtain the functions a^(q) we use equation (9). However, in some cases which we prefer to 
use the Hamiltonian formaUsm of the theory, the Hamiltonian counterpart of this equation is better 
suited in finding these coefficients which, equivalently, can be written as 

{Pa, H} + ^ {q'',h} Pa = 0. (11) 

In general, the expression above gives a quadratic polynomial in terms of momenta with coefficients 
being partial derivatives of with respect to the configuration variables q. Thus, the expression 
is identically equal to zero if and only if these coefficients are zero, leading to a system of partial 
differential equations for Q;"^(q). An important ingredient in any model theory related to Noether 
symmetry is the choice of the cyclic variables related to this symmetry. Since such a variable gives 
a constant of motion, it should have a vanishing commutator with the Hamiltonian and therefore 
from a quantum mechanical point of view they should have simultaneous eigcnfunctions. This means 
that to describe the quantum structure for a Noether symmetric model, one should take into account 
the commutative algebra between the Hamiltonian and the conserved quantities (10). If Q is itself a 
cyclic variable the quantum counterpart of the theory can be described by the following equations 

' H*(q) = 0, 

(12) 

^ Q^r(q) = Q^(q), 

where Q is the operator form of (10) and Q is its eigenvalue. If Q is not a cyclic variable this procedure 
does not work. In this case we seek a change of variables in the form of a point transformation 
q = (q^jQ^,---) — ?• u = (ti^,n^,...) on the vector field (8) such that in terms of the new variables u, 
the Lagrangian includes one (or more) cyclic variable. A general discussion of this issue can be found 
in [11]. Under such a point transformation it is easy to show that the vector field (8) takes the form 
[13] 

it is easy to show that if X is a Noether symmetry, X has also this property, that is, XjC. = 
XjC. = 0. Thus, if we demand Xu^ = 1 for some G u and Xu^ = for any j ^ i, we get 

X = ^^X£ = ^ = 0. (14) 

This means that is a cyclic variable and the dynamics can be reduced. On the other hand, the 
constant of motion Q becomes 

Q = a-ic)PA = a-i^)^ = a-i^) ^-^^ + j . (15) 

f) B Pi' B ft B 

Since q ^ u is a point transformation, we have 1^ = and = f^. Therefore, 

^ = " = ^-'di^B = 5^ = (16) 

where tTj is the momentum conjugate to u^. Thus, as expected the constant of motion which cor- 
responds to the Noether symmetry is nothing but the momentum conjugated to the cyclic variable. 
Now, it is obvious that [iri^H] = and the quantum description of the system under consideration 
can be given by the following equations 

H*(q) = 0, 

(17) 

7fi*(q) = C7i^'(q), 
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where (Tj is the eigenvalue of vfj. 

In the next section we will apply the above formalism to a given cosmological model as a dynamical 
system. In cosmological systems, since the scale factors, matter fields and their conjugate momenta 
play the role of dynamical variables, introduction of Noether symmetry by adopting the approach 
discussed above is particularly relevant. 



3 Scalar field cosmology 



In this section wc consider a FRW cosmology with a scalar field with which the action of the model 
is augmented. In a quasi-spherical polar coordinate the geometry of such a space-time is described 
by the metric 



-N\t)dt^ + a\t) 



1 — kr'^ 



+ r 



a2 , -2 
^ + sm 



(18) 



where N{t) is the lapse function, a(t) the scale factor and A; = 0, ±1 is the curvature index. Since 
our goal is to study the models which exhibit Noether symmetry, we do not include any matter 
contribution in the action. Let us start from the action (we work in units where c = % = IGttG = 1) 



= j d^X^/^ 



(19) 



where g is the determinant of the metric tensor g^y^ R is the Ricci scalar corresponding to g^y^ y{(f) is 
the potential function for the scalar field <j){t) and the parameter e = ±1 corresponds to the ordinary 
scalar field (where e = +1) or phantom scalar field (where e = —1). By substituting (18) into (19) 
and integrating over spatial dimensions, we are led to a point-like Lagrangian in the minisuperspace 
{a, ^} as 

C = -Mh^ + ]^ea^4p' + 2,ka- a^V{cl)), (20) 



in which we have set N = 1 so that the time parameter t is the usual cosmic time. Now, it is easy to 
see that this minisuperspace has the following minisupcr metric 



QAsdq dq 



B 



-6ada^ + ea^dcl)^ . 



(21) 



This two dimensional minisuperspace is the space spanned with the dynamical variables of the model 
and indeed is the space in which we are looking for a Noether symmetry. In the following we will apply 
the Noether symmetry approach to the minisuperspace which is represented by a curved manifold 
with a minisuper metric given by (21). 



3.1 Ordiiiciry scalar field 

In this case (e = +1) the minisuperspace has a Lorentzian metric. However, a simple calculation shows 
that its Ricci scalar is zero and therefore this space is a two dimensional Lorentzian flat space. This 
means that there is a suitable change of variables that can transform the metric into a Minkowskian 
one. To do this, consider the following change of variables q'^ = (a, (j)) — >■ = {x, y) [14] 

X = -a^/^ cosh(a;^), y = -a^/^ sinh(w^), (22) 

LJ CO 

where uo = s/G/A. In terms of these new variables, Lagrangian (20) takes the form 

^=l{i'- - - V') y (f ) + 3fcc.^/^ (x' - y') , (23) 

with the corresponding Hamiltonian becoming 

W = \ {Pl - P^) + - y') V (I) - ?>ku'" {x' - y') , (24) 



5 



where and Py are the momenta conjugate to x and y respectively. Thus, it is easy to see that in 
the minisuperspace constructed by Q"^ = {x,y), the metric is Minkowskian and represented by 



GAEdQ^'dQ^ = -dy' + dx\ 



(25) 



Now, we have a set of variables (x^y) endowing the minisuperspace with a Minkowskian metric 
and we are going to find the potential function V{y/x) such that the dynamical system described 
with the Lagrangian (23) exhibits the issue of the Noether symmetry. For this purpose, following the 
steps after equation (8) we define the generator of the desired Noether symmetry as the vector field 



d n d . d ■ d 
X = a—+P— + a—+ /3— , 
ox oy ox oy 



(26) 



on the tangent space of the corresponding configuration space such that the Lie derivative of the 
Lagrangian with respect to this vector field vanishes, that is, LxC = 0. In (26) a and /3 are some 
functions of x and y which in order to obtain them we use the Noether symmetry condition. For 
Lagrangian (23) this condition results in 



a 



u:^xV{z) + ^Lo^^ (x2 - y^) V'iz) + 2kco^/-'x {x^ - y^) 



-2/3' 



+ 



u:''yV{z) 
da . da 



1 2I 

-UJ — 

2 X 



dx ^ dy 



x^-y^) V'{z)-2ku''/^y [x'' - y^) 

d(3 



-2/3 



+ 



_■ 
dx dy^ 



y = o, 

where z = y/x and V'{z) = dV/dz. Now, we are led to the following system of equations 



(27) 



L0^Viz)i(3y - ax) + ^wV(z) (x^ - y^) (^a^ - /3^) + 2kuj^'^ - y^) '^'^ {ax - ^y) = 0, (28) 

dx dy dy dx 

Equations (29) admit the independent solutions as 

(a,^) = (l,0), (a,^) = (0,l), {a,P) = {y,x), (30) 

which upon substitution into relation (28) yields the potential function. In what follows we restrict 
ourselves to the flat case k = and consider three cases separately. 

3.1.1 The case: (a,/?) = (1,0) 

This solution is related to the Noether symmetry generator X = and Noether conserved charge 
Q = Px- In this case from equation (28) we obtain 

V('-)-^- (31) 

\x J — 

In terms of the old variables (a, ^) this potential is ^^(0) ~ sinh^ (a;(?!)) . With the choose of this 
potential the Hamiltonian of the model becomes 

n = l{p'x-P^)-l^V. (32) 

This Hamiltonian is nothing but the difference of a free particle and a harmonic oscillator Hamilto- 
nians. Classically, this means that the motion in y-direction is bounded and quantum mechanically 
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it becomes natural to require the corresponding wave function to be normalizable. Now, the classical 
and quantum solutions of the model described by Hamiltonian (32) can be easily obtained. The 
classical dynamics is governed by the Hamiltonian equations, that is 



X = {x,n} = P:c, = {P^^n} = 0, 

[ y = {y,n} = -Py, Py = {Py,n} 



These equations can be immediately integrated to yield 

X{t) = Poxt + Xo, Px = Pqx, 

and 



y(^t) = — sin (cot + 60), Py{t) 



-Pqx cos {ut + 5q) , 



(33) 



(34) 



(35) 



where Pq^^ xq and 5q are integration constants. Since in the quantum version of the model we are 
interested in constructing wave packets from the WDW equation, we would like to obtain a classical 
trajectory in configuration space {x,y), where the classical time t is eliminated. This is because no 
such parameter exists in the WDW equation. It is easy to see that the classical solutions (34) and 
(35) may be displayed as the following trajectories 



sm 



UJ 



-—{x-xo) + 5q 



(36) 



Going back to the old variables (a, (f)) we get the cosmological solutions from (22) as 

1/3 



a{t) = UJ 



= ,,,2/3 



p2 

{Poxt + xof - sin2(a;t + do) 



and 



(p{t) = — Arctanh 

UJ 



^sin {ujt + 60) 

Poxt + Xq 



(37) 



(38) 



We now focus attention on the study of the quantization of the model described above. This can be 
achieved via canonical quantization procedure which leads the WDW equation. For the Hamiltonian 
(32) this equation reads 

/ 52 a2 \ 

-^-u;V]^{x,y) = 0. (39) 



dx'^ dy^ 



We separate the variables in this equation as "^{x, y) = e'''^^y{y) leading to + {v^ — uj^y^)y{y) = 0, 
where u is a separation constant. The physical acceptable solutions to this equation can be found 
in terms of the Hermitc polynomials Hn{z), if u'^ = (2n + l)uj = I'^w, for some integer n. Thus, the 
eigenfunctions of the WDW equation can be written as 



*„(x,y) = e-"^e-^^'/fn(y). 



(40) 



where (X, Y) = ^/uj(x, y). There still remains the question of the boundary conditions on the solutions 
to the WDW equation. Note that the minisuperspace of the above model is a two dimensional manifold 
< a < 00, —00 < (j) < +00 . According to [15], its nonsingular boundary is the line a = with 
< 00, while at the singular boundary, at least one of the two variables is infinite. In terms of 
the variables x and y, introduced in (22), the minisuperspace is recovered by x > 0, x > and 
the nonsingular boundary may be represented by x = y = 0. Since the minisuperspace variables 
arc restricted to the above mentioned domain, the minisuperspace quantization deals only with wave 
functions defined on this region. Therefore, to construct the quantum version of the model one should 
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take into account this issue. This is because that in such cases one usually has to impose boundary 
conditions on the allowed wave functions otherwise the relevant operators, specially the Hamiltonian, 
will not be self-adjoint. The condition for the Hamiltonian operator % associated with the classical 
Hamiltonian function (24) to be self-adjoint is {il^i,'Hip2) = (^■0i,V'2) or 

/ Vi{x,y)iiip2{x,y)dxdy = I ip2{.x,y)ii'il)l{x,y)dxdy, (41) 

where the integrals should be taken over the domain where the minisupcrspacc variables are defined on 
which. Following the calculations in [16] and dealing only with square integrablc wave functions, this 
condition yields a vanishing wave function at nonsingular boundary of the minisuperspace. Hence, we 
impose the boundary condition on the solutions (40) such that at the nonsingular boundary (at a = 
and < oo) the wave function vanishes. This makes the Hamiltonian hcrmitian and self-adjoint 
and can avoid the singularities of the classical theory, i.e. there is zero probability for observing a 
singularity corresponding to a = 0. Therefore, we require 

^(a = 0, 0) = ^ ^'(X = 0, y = 0) = 0, (42) 

which yields 

Hn{0) = 0^n = odd. (43) 

In general, one of the most important features in quantum cosmology is the recovery of classical 
cosmology from the corresponding quantum model or, in other words, how can the WDW wave 
functions predict a classical universe. In this approach, one usually constructs a coherent wave 
packet with good asymptotic behavior in the minisupcrspacc, peaking in the vicinity of the classical 
trajectory [17]. Therefore, we may now write the general solutions to the WDW equation as a 
superposition of the eigenfunctions, that is 

oo 

^{X,Y) = ^ c„e-"^e-^^'iJ2n+i(l^), (44) 

n=0 

where Cn are suitable weight factors to construct the wave packets and to make a convergent su- 
perposition, we may choose them such that the summands with smaller n have more important 
contribution to the above superposition. Therefore, to achieve an analytical expression for the wave 
function we assume that the above superposition is taken over such values of n for which one can use 
the approximation i/„ ~ 2n -|- 1, so that 

oo 

*(X, Y) = J2 cne^('"+^)^e-i^'//2n+i(l^). (45) 

n=0 

By using the equality 

E 7^^rXW^2n+i(^) = e-^ sinh272;, (46) 

we can evaluate the sum over n in (45) and simple analytical expression for this sum is found if we 
choose 7 to be e*^ and c„ = J\f/{2n + 1)!, which results in^ 

*(X, Y) = Afe~^^^ exp(-e2^^) sinh (2Ye'^) , (47) 

where M is a numerical factor. Before going any further, some remarks are in order. An important 
ingredient in any model theory related to the quantization of a cosmological setting is the choice 



^One may have some doubts on this final form for the wave function and the following results due to the assumption 
Vn ^ 2n + 1 which seems to be irrelevant especially in the case of a discrete spectrum for n. To overcome this problem, 
we have made a numerical study of the behavior of based on equation (44) with c„ = l/(2n-|- 1)! and Un = {2n+l)^ 
and have verified that the general patterns of the resulting wave packets follow the behavior shown in figure 1 with a 
very good approximation. 
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Figure 1: The figures show |\E'(X, y)p, the square of the wave function and its corresponding contour plot. Up: the 
figures are plotted for negative values of y and bottom: the same figures are plotted for positive values of y. 

of a probability measure in order to make predictions from a given solution to the WDW equation. 
Since the WDW equation is a second order Klaein-Gordon type equation in quantum field theory, 
the most widely used probability density for it results the conserved current j = |(^'*V^' — ^'V^'*). 
But this immediately presents a problem, because the probability density given by this expression 
is not positive definite. For this reason and some other difficulties related to this interpretation, an 
alternative used approach is to introduce l^'p as the probability density. Although this definition is 
necessarily positive definite and specially works very well for the minisuperspace models, it suffers 
from some problems. The most important of such problems arise when we are dealing with the 
non-normalizable. Despite of the existence of the above mentioned problems, here we shall look for 
the peaks in the wave function in the sense that if a peak is sufficiently strong in some region, the 
probability of the universe being in this region is large. In figure 1 we have plotted the square of 
the wave function (47) and its corresponding contour plot. It is seen that the peaks of the wave 
function follow a periodic pattern. In comparison with the classical path (36), this means that a good 
correlation exists between the quantum patterns shown in this figure and the classical trajectory (36) 
in configuration space {x,y). 

3.1.2 The case: (a, ^) = (0,1) 

This case relates the Noether vector ^ and its corresponding conserved charge is Q = Py. From (28) 
the potential function reads as 




where in terms of the old variables (a, cf)) can be termed as V{4>) ~ cosh^ Since the Hamiltonian 

of the model, in this case is 

^='^{P'-Py)+r'^^ (49) 

we see that there is no major difference between this case the the previous one except that the 
variables x and y change their role with each other. Therefore, following the steps (32)-(47), we are 
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led to the classical solutions 

x{t) 



Poy_ 



sin {u)t + r/o) , Px = Poy cos {ut + r/o) , 



yit) = -Poyt + yo, Py{t) = Poy, 
and the corresponding classical trajectories 



Pay . 

X = — - sm 

U! 



00 



-f^{yo - y) + m 



Oy 



(50) 



(51) 



(52) 



where in terms of the old variables (a, 0) can be viewed as 

a(t)=a;2/3 



P? 



^ siii^{ut + r]o) - {-Poyt + yof 



1/3 



and 



-Ppyt + yo 



.^sin (oot + m). 



(f>{t) = — Arctanh 

CO 

Also, the wave function of the quantum counterpart of the model can be written as 

^{X, Y) = J\fe-^^^ exp(-e2^^) sinh (2Xe^^) . 



(53) 



(54) 



(55) 



The discussions on the comparison between quantum cosmological solution and its classical version 
are the same as previous model, namely the (a, /3) = (1, 0) model. Similar discussion as above would 
be applicable to this case as well. 

3.1.3 The case: (a, /3) = (y, 

In this case the Noether symmetry is generated by the following vector field 



d d d d 
ox ay ox oy 



(56) 



which corresponds to the Noether conserved charged Q = yP^ + xPy. The condition (28), in this 
case, demands that the potential function be a constant function and as such, we choose it to be 1. 
Therefore, the Hamiltonian takes the form 



(57) 



which describes an isotropic oscillator-ghost-oscillator system. The classical equations of motion are 
given by 



x = {x,n} = Px, Px = {Px,n} = -uj^x, 

y = {y,H} = -Py, Py = {Py,H} = co^y. 



(58) 



Choosing the integration constants such that the solutions satisfy the zero energy condition % = 
the solutions are obtained as 



x{t) = A sin {ujt + rjo) , Px{t) = Aco cos {oot + rjo) , 
y{t) = £Asm {cot + Sq) Py{t) = -£Alo cos {cot + do) 



(59) 
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where A, 6q and r]o are integration constants and i = ±1. From the above equations, we see that the 
classical trajectories obey the relation 

y'^ + - 2£xy cos(5o - Vo) - sin2((5o - r/o) = 0. (60) 

This equation describes ellipses which their major axes make angle 7r/4 with the positive/negative 
X-axis according to the choices ±1 for I. Also, the eccentricity and the size of each trajectory 
are determined by {5q — r/o) and A respectively. Now, using the transformation (22) the classical 
cosmological behavior can be obtained as 

a{t) = A2/3a;2/3 sni^/^{5Q - r/o) sin^/^ {2ojt + 5q + rjo) , (61) 

and 

?sin [cot + 5q) 



6(t) = — Arctanh 



(62) 

sm [iot + ryo) J 

At this step, as the previous subsections, we deal with the quantization of the model. The WDW 
equation corresponding to the Hamiltonian (57) reads 

(-^ + ^ + - ^"y') *(^' y) = 0- (63) 

This equation is a quantum isotropic oscillator-ghost-oscillator system with zero energy. Therefore, 
its solutions belong to a subspace of the Hilbcrt space spanned by separable cigcnfunctions of a 
two-dimensional isotropic simple harmonic oscillator Hamiltonian. In [18] the models of coupled 
harmonic oscillators have been studied in other contexts relevant for quantum cosmology. Separating 
the eigenfunctions of (63) in the form 

^m,n, {X, Y) = {X)yn, (Y), (64) 

yields 

Xm{X)=e--^'''HnAX), :y„,(y) = e-^^'if„,(y), (65) 

subject to the restriction n\ = n2 = n. In (65), Hn{z) are the Hermite polynomials and {X,Y) = 
\foj(x^ y) as before. Applying again the boundary condition (42), we are led to the following general 
solution to the WDW equation 

oo oo 

^(X,Y) =J2Y. c„,c„,e-5(^'+^')i72„^+i(X)iJ2n2+i(^)- (66) 

ni=0 n2=0 

If we choose the coefficients 0^ and c„2 to be Cm = ^'^^^'^^ / {2ni + 1)! and c^j = C^"^'''^/(2n2 + 1)! 
where ^ and are arbitrary complex constants, the wave function takes the form 

*(X, Y) = Are-5(^'+'^')e-(«'+^') sinh (2^X) sinh {2CY) . (67) 

In figure 2, we have shown the qualitative behavior of this wave function and its contour plot. To 
realize the correlation between these quantum patterns and the classical trajectories represented by 
(60), note that in the minisuperspace formulation, the cosmic evolution of the universe is modeled 
with the motion of a point particle in a space with minisuperspace coordinates. In motion on an 
ellipse, the particle (universe) spend most of its time in the region near its apogees. This means that 
the particle will be found around the apogees of its trajectory with the maximum probability. This 
is just what the figure 2 is showing. As it is seen from this figure, the wave function has two pair 
peaks, each pair may be interpreted as two apogees of the classical trajectory. In classical model, 
the transition from one point on the configuration space to any other points will be described by a 
continuous motion on the classical trajectory. The quantum description of such a transition, on the 
other hand, may be explained by a tunneling procedure. This means that there are different possible 
states from which our present universe could have evolved and tunneled in the past, from one state 
to another. Again, we see that our quantization leads us to a model free of the classical singularity 
and a good correlation with its classical counterpart. 
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Fi gure 2: The figures show the square of the wave function (67), its corresponding contour plot and the 

classical trajectory (60), (solid line for £ — 1 and dashed line for i = —1). 



3.2 Phantom scalar field 

In this case (e = —1), the minisuperspace manifold would have an Euclidean signature (see (21)) with 
zero Ricci scalar. To write the metric in its canonical form, we apply the transformation 

X = — o^/^ cos(a;(/)), y = —a^^"^ sm{uj<f>). (68) 

U} UJ 

It is easy to see that in terms of these new variables the metric takes the form of a 2-dimensional 
Euclidean one represented by 

GABdQ^dQ^ = -dx'^ - dy^. (69) 
Also, the phantom counterpart of Lagrangian (20) becomes (we consider k = 0) 

C = \{-i'-i^)-l.'(.'^y')v(i). (70) 

with the corresponding Hamiltonian becoming 

Again, considering a Noether symmetry generated by a vector field like (26), the condition LxC = 
will lead us to the following equations for a and /3 

-uj^V{z){(3y + ax) + \uj^V\z) (x^ + y^) [a^ - = 0, (72) 

^ = ^ = 0, ^ + ^ = 0, (73) 
dx dy ^ dy dx ' 

where z = y/x and V'{z) = dV/dz as before. The independent solutions of the system (73) my be 
obtained as 

(a,/3) = (1,0), (a,/3) = (0,1), (a, /3) = {y, -x). (74) 

Now, following the same procedure as in the previous subsection, we consider the above tree case 
separately to study their corresponding classical and quantum cosmology. 

3.2.1 The case: (a, /3) = (1,0) 

This case corresponds to the vector field X = which results the conserved quantity Q = Px- The 
potential function can be obtained from (72) as 
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where in terms of the cosmic variables (a, (6) may be written as V{<^ ~ sin^ (a;(?!)) . This potential 
function transforms the Hamiltonian (71) into the form 



2 2 

w y , 



2 V ^ yj ■ 2 

from which one gets the following classical equations of motion 

' x = {x,n} = -p^, = {p^, n} = 0, 



(76) 



with solutions 
and 



[y = {y,n} = -Py, Py = {Py,H} = -u^y, 



X{t) = -Poxt + Xq, Px = Pox, 



y(t) = cosh {ojt + 6o) , Py{t) = -Pqx sinh {ujt + Sq) , 



(77) 



(78) 



(79) 



where for the sake of the simplicity, we use the same notation as in the previous section for the 
integration constants. These solutions obey the classical trajectories in the configuration space {x,y) 



y = cosh ( — (zo - x) + 6o 



(80) 



Finally, going back to the original variables {a,(l)), the corresponding classical cosmology will be 
obtained as 



and 



ait) = a;2/3 



p2 

{-Poxt + Xof + ^ COSh^ {ujt + So) 



1/3 



1 



Arctan 



CO 



^ cosh (ujt + So) 

-Poxt + Xo 



(81) 



(82) 



We see that unlike the case of the ordinary scalar field, here all the motions are unbounded, so that 
it becomes natural to expect that the wave function describing the quantum version of the model to 
be non-normalizable. Now, by canonical quantization the WDW equation is given by 



&2+^2+-v)nx,y) = o. 



(83) 



The eigenfunctions of this equation can be found by the method of separation of variables, by means 
of which we obtain 



*,(X,r) = e±-^P_i 2 {{±l + i)Y), 

2=^* 2 



where is the separation constant and I? is the parabolic cylinder function. As before the general 
solutions may be constructed by a superposition of these eigenfunctions. Since non-normalizable states 
in this case belong to the continuous spectrum, the general wave function should be constructed by 
a continuous superposition in the form 



*(x,r) 



+ 00 



C{iy)^^iX,Y)du. 



(85) 



However, to achieve a more clear picture of the wave function by the numerical methods which would 
be comparable with the results in the previous section, we restrict ourselves to a discrete superposition 
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Figure 3: From left to right: the square of the phantom wave function (86), in which we have taken the 

combinations from six terms (taking more terms would only have minor effects on the shape of the wave function) with 
c„ = l/(n!2"), its corresponding contour plot and the classical trajectory (80). 

of the above eigenfunctions. Thus, by choosing the separation constant as = (2n — i)uj = v^uj and 
connecting the parabohc cyhnder function with the Hermite functions [19], the general solution to 
the WDW equation takes the form 

oo / 1 \ 

*(X,y) = ^ c„e±-"^e-5^^'//,„ —(1 -i)Y\. (86) 

We have summarized the above results in figure 3. As this figure shows, like in the case of the 
ordinary scalar field, the evolution of the universe based on the classical cosmology represents a late 
time expansion. This is because the particle (universe) moves on one of the branches of its classical 
trajectory forever without any turning points. A remarkable point about this motion is that in spite 
of the ordinary scalar field case, it does not come from a big-bang singularity in which the scale factor 
goes to zero, but instead tends to another kind of singularity in which the scale factor diverges. In 
the models where the phantom fields are considered by a perfect fluid with the equation of state 
parameter less than —1, this kind of singularity will be achieved in a finite-time and is called the 
big-rip singularity. On the other hand a glance at the quantum patterns shows that although a good 
correlation exists between them and the classical loci in the configuration space, but the quantum 
effects dominate in the region of the classical singularity, i.e., at the large values of scale factor. In this 
regime the quantum solutions fall from an expansion phase to a contraction era and this phenomenon 
will be repeated cyclic. By such a behavior we see that the quantum effects show their important 
role for the large values of the scale factor and thus the evolution of the scale factor towards a big-rip 
like singularity will be avoided [20]. However, a remark about the above analysis is that it happens 
in some internal time parameter such as x and in order to see whether it is physically meaningful the 
results should translate in terms of the cosmic time t. The answer to this question in general may 
be not an easy task. This is because one usually uses non-standard parametrization of the metric 
in order to simplify the calculations and have a manageable Lagrangian. Therefore, returning the 
results into the proper time gauge is not always an integrable process. By the way, in our simple 
model at hand, since (78) represents a monotonic relation between x and t, we may argue that the 
above conclusions are also valid in terms of the cosmic time gauge t as well. 

3.2.2 The case: (a, /3) = (0, 1) 

The situation in this case is alike the one that was occurred when we were discussing about the 
ordinary scalar field. The potential function resulting from this solution is y(y/x) = x^ jix? -|-y^), or 
y((/)) ~ cos^(w(^) in terms of the variables (a, 0). Again, all of the results coming from this potential 
function can be obtained from their counterparts in the previous subsection with exchange the role 
of the variables x and y with each other and because of this we do not repeat the calculations ones 
again. 
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3.2.3 The case: (a,/3) = 

This solution generates the Noether symmetry with the vector field 

d d d d 
ox oy ox Oy 

which gives the conserved charge Q = yP^ — xPy. Substitution of this solution into the equation (72) 
we get a constant value for the potential function and as before we take it to be equal to 1. Therefore, 
the corresponding dynamical system can be described by the Hamiltonian 



2 V ^ yj 2 

Following the same steps as in the above sections, we arrive the classical equations of motion as 



(88) 



X = {x^n} = -p^, p^ = {p^,n} = -oj^x, 
[y = {y,n} = -Py, Py = {Py, n} = -u^y, 

which admit the following integrals 



(89) 



x{t) = A cosh {cot + r)o) , Px{t) = —Auj sinh {ojt + rjo) , 

y{t) =£Asmh{ujt + 6o) Py{t) = -Mw cosh (wt + (5o) , 



(90) 



where as before we take A, 5q and rjo as integration constants and £ = ±1. Now one can eliminate 
the time parameter t from these solutions to get the classical trajectories in the configuration space 
(x, y) as 

-x^ - 2£xy sinh((5o - r/o) + A^ cosh^((5o - m) = 0, (91) 

which represents a hyperbola so that the the particle (universe) has an unbounded motion in the 
x — y plane. Now, going back to the original variables (a, (f)) with the help of the transformation (68), 
the dynamics of the system reads as 



and 



a{t) = J^l^u?!^ coshV3(5o _ j^^) coshes (2a;t + 5o + r/o) , 

"£sinh {oot + (5o) 



(?!)(t) = — Arctan 



(92) 



(93) 



cosh [lot + ?7o) 

To proceed, let us to deal with the quantum cosmology associated with the model described by the 
Hamiltonian (88). The corresponding WDW equation is 



' q2 q2 \ 



(94) 



With the same techniques as in the previous sections and also with same notations one may represent 
the general solutions to the above equation as 



^(x,y) = ^c„e-^ 

n=0 



(-^(1 - i)X^ //-i-in (^(1 -i)Y\ + 

H^n (^(1 - i}Y'^ i^-l-m (^(1 " i)X 



(95) 
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Figure 4: From left to right: the square of the phantom wave function (95), in which we have taken the 

combinations from six terms (taking more terms would only have minor effects on the shape of the wave function) with 
c„ = l/(n!2"), its corresponding contour plot and the classical trajectory (91) with £ = 1. 

Figure 4 shows the quahtative behavior of the wave function and classical trajectory for typical values 
of the parameters. As we have mentioned above, the motion in configuration space is an unbounded 
motion along one of the branches of a hyperbola. A point (universe) with such a motion spend lots 
of its time far from the apex of its path. This is exactly what we realize from the quantum wave 
function in figure 4. These patterns show that the peaks of the wave function follow the classical 
trajectory with a good degree of accuracy. 

4 Conclusions 

In this paper we have studied the scalar field classical and quantum cosmology in a Noether symmetry 
point of view. The minisuperspace of such a model is constructed of a two dimensional manifold with 
vanishing Ricci scalar. Therefore, it is possible to find a coordinate transformation which cast the 
corresponding metric to a Minkowskian or Euclidean one according to the choices of an ordinary or 
phantom model for the scalar field. Here, after a brief review of the issue of the Noether symmetry in 
a dynamical system we have applied the formalism into a FRW cosmological model with a scalar field 
as its source. The phase space was then constructed by taking the scale factor a{t) and scalar field 
(p{t) as the independent dynamical variables. The Lagrangian of the model in the configuration space 
spanned by {a, (/>} is so constructed such that its variation with respect to these dynamical variables 
yields the appropriate field equations. As is mentioned above, we introduced a new set of variables 
{x, y) in terms of which the minisuperspace takes the form of a Minkowskian space, in the case of a 
ordinary scalar field, and a Euclidean space in the case where the scalar field is of phantom type. Also, 
the dynamics of the (x, y) system is described by a Lagrangian which is a diagonal quadratic form 
with constant coefficients. The existence of Noether symmetry implies that the Lie derivative of this 
Lagrangian with respect to the infinitesimal generator of the desired symmetry vanishes. By applying 
this condition to the Lagrangian of the model, we obtained the explicit form of the corresponding 
potential function. In the both cases of ordinary or phantom scalar field we found three distinct 
solutions for the potential function each of which have their own classical and quantum dynamics. 
In more details , in the case where the scalar is an usual one the three forms of Noether symmetric 
potentials are found as V{y/x) = y'^/{y'^ — x^), V{y/x) = — y^) and V = cons. = 1. In the 

first two cases one of the variables has a cyclic dynamics while the other behaves linearly as time 
progresses. Also, for a constant potential function, we saw that the system of the two variables (x, y) 
is a two dimensional oscillator-ghost-oscillator system in which both of the variables oscillate with the 
same frequency. The classical trajectories of these models are also obtained and it is seen that they 
all represent a bounded motion in the configuration space. We have shown that all of these classical 
solutions have a singularity of the big-bang or big-crunch type and so to pass this issue we have dealt 
with the quantization of the model. As for the quantum version of these models, we obtained exact 
solutions of the WDW equation. The wave function of the corresponding universe consists of some 
branches where each may be interpreted as part of the classical trajectory. We saw that since the 
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peaks of the wave function follow the classical trajectory, there seems to be good correlations between 
the corresponding classical and quantum cosmology. The same study for a phantom scalar field yields 
the potentials V{y/x) = /{x^ + y^), V{y/x) = x'^/{x^ + y^) and V = 1. For these potentials, we 
solved the classical equations of motion exactly and showed that at least one of the variables behaves 
hyperbolically. This causes the classical trajectory in x — y plane to be an unbounded path in which 
the moving particle (universe) tends to regions where the scale factor is large. The resulting quantum 
cosmology and the corresponding WDW equation in the phantom framework were also studied and 
analytical expressions for the wave functions of the universe were presented with good correlations 
with the classical trajectories. 

Finally, we would like to emphasize that for an ordinary scalar field the quantum effects dominate 
in the region of the classical big-bang singularity, i.e., at the small values of scale factor. At the 
big-bang the quantum solutions bounce from a contraction phase to an expansion era. On the other 
hand for a phantom field the quantum effects dominate in the region of the classical large scale factor 
so that in this region quantum solutions fall from an expansion phase to a contraction era. 
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